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ABSTRACT 

Context. The findings of more than 350 extrasolar planets, most of them nontransiting Hot Jupiters, have revealed correlations between 
the metallicity of the main-sequence (MS) host stars and planetary incidence. This connection can be used to calculate the planet 
formation probability around other stars, not yet known to have planetary companions. Numerous wide-field surveys have recently 
been initiated, aiming at the transit detection of extrasolar planets in front of their host stars. Depending on instrumental properties 
and the planetary distribution probability, the promising transit locations on the celestial plane will differ among these surveys. 
Aims. We want to locate the promising spots for transit surveys on the celestial plane and strive for absolute values of the expected 
number of transits in general. Our study will also clarify the impact of instrumental properties such as pixel size, field of view (FOV), 
and magnitude range on the detection probability. 

Methods. We used data of the Tycho catalog for x 1 million objects to locate all the stars with 0'" < my < 11.5'" on the 
celestial plane. We took several empirical relations between the parameters listed in the Tycho catalog, such as distance to Earth, 
my, and (6 - V), and those parameters needed to account for the probability of a star to host an observable, transiting exoplanet. 
The empirical relations between stellar metallicity and planet occurrence combined with geometrical considerations were used to 
yield transit probabilities for the MS stars in the Tycho catalog. Magnitude variations in the FOV were simulated to test whether this 
fluctuations would be detected by BEST, XO, SuperWASP and HATNet. 

Results. We present a sky map of the expected number of Hot Jupiter transit events on the basis of the Tycho catalog. Conditioned by 
the accumulation of stars towards the galactic plane, the zone of the highest number of transits follows the same trace, interrupted by 
spots of very low and high expectation values. The comparison between the considered transit surveys yields significantly differing 
maps of the expected transit detections. While BEST provides an unpromising map, those for XO, SuperWASP, and HATNet show 
FsOV with up to 10 and more expected detections. The sky-integrated magnitude distribution predicts 20 Hot Jupiter transits with 
orbital periods between 1.5 d and 50 d and my < 8™, of which two are currently known. In total, we expect 3412 Hot Jupiter transits 
to occur in front of MS stars within the given magnitude range. The most promising observing site on Earth is at latitude = - 1. 

Key words. Stars: planetary systems - Occultations - Galaxy: solar neighborhood - Galaxy: abundances - Instrumentation: miscel- 
laneous - Methods: observational 

1. Introduction prototype, was unknown and previous considerations had been 

H . 11 . based on syster ns similar to the so lar system. Using geometrical 

.C^ , A short essay by Otto Struve (|Struve||l952|) provided the first considerations. iRosenblatj [19TW\ found that the main contri- 

published proposal of transit events as a means of exoplane- bution to the transit probability of a solar system planet would 

tary detection and exploration. Calculations for transit detec- come from the inner rocky planets. However, the transits of these 

tion probabilities (R osenblatt 197 1 ; Borucki & Summers 1984; relatively tiny objects remain undetectable around other stars as 
iPepper & Gaudil2006l) and for the expected properties of the dis- 

covered plaiiets have been done subsequently by many others ^^^^ transit of exoplanet was finally detected 

(iGiUon et alj|2005t [Fressin etalj|2007|: |Beatty & Gaudi||2008|). around the sun-like star HD209458 (Charbonneau et al. 200(11: 



around the sun- like star HD209458 (Charbonneau et al. 200Q|; 
Until the end of the 1990s, when the sample of known exoplan- [pueloz et al.l I2OOOI) . Thanks to the increasing number of ex- 
ets had grown to more than two dozen (j Castellano etal. , 2000D , opi^net search programs, such as the ground- based Optical 
the family of so-called 'Hot Jupiters', with 51 Pegasi as their Gravitational Lensing Experiment (OGLE) (lUdalski et al 



*Sky maps (Figures 1 and 3) can be downloaded in electronic form '2002 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 



1992), the Hu ngarian Automated Telescope (HAT) ( IBakos et al 



|2004 . the Super Wide Angle Search for Planets 



via,http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-HA/| 'A correction to his Eq. (2) is given in lBorucki & SummersI ( Il984h . 
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(SuperWASP) (IStreet et al.ll2003h. the Berlin Exoplanet Search 
Teles cope (BEST) (iRauer et al.ll20()4l) . XO dMcCullough et alj 



2005), the Transatlantic Exoplanet Survey (TrES) ( Alonso et aL 
2007h . and the Tautenburg Exop lanet Search Telescope 



(TEST) jEigmiiller & Eisloffelll2009l) and the space-based mis- 
sions 'Convection, Rotation & P lanetary Transits' (CoRoT) 
('Baglin et al.l l2002h and Kepler (IChristensen-Dalsgaard et all 
[2007), the number of exoplanet transits has grown to 62 until 
September 1^' 2OOS0 and will grow drastically within the next 
years. These transiting planets have very short periods, typically 
< 10 d, and very small semimajor axes of usually < 0.1 AU, 
whic h is a selection effe ct based on geometry and Kepler's third 
law (iKepler et al.lll619T) . Transiting planets with longer periods 
present more of a challenge, since their occultations are less 
likely in terms of geometrical considerations and they occur less 
frequently. 

Usually, authors of studies on the expected yield of tran- 
sit surveys generate a Active stellar distribution based on stel- 
lar population models. iFressin et al.l ( l2007b use a Monte-Carlo 
procedure to synthesize a Active stellar field for OGLE based 
on star counts from Gould et al. (2006^, a stellar metallicity dis- 



tribution from iNordstrom et al. (|2004|). and a synthetic struc 



ture and evolution model of Robin et aP (l2003h . The metallicity 
correlation, however, turned out to underes timate the true stel- 
lar metallicity by about 0. Idex, as found byl^antos et al." ("20041) 
andjFischer & Valenti (2005,). In their latest study, Fress in et aP 



(l2009l) first g enerate a ste l lar po pulation based on the Besan§on 
catalog from |Ro bin et al. ( 2003h and st atistics for multiple sys- 
tems from .Duquennov & Mavoi J 199l|) to ap plv then the metal- 
licity distri bution from Santos et al. (2004) and issues of de- 
tectabihty (iPont et alj [2006) . Beatty & Gaudi ( 2008i) r elv on a 



Galactic structure model bv iBahca 



& Soneiral (1198^ 

function as suggested bv lReid et al.l (l2002i) based on Hipparcos 
data, and a model for interstellar extinction to estimate the over- 
all output of the current transit surveys TrES, XO, and Kepler 
In their paper on the number of expected planetar y transits to b e 
detected by the upcoming Pan-STARRS survey ( lKaiseiil2004l) . 
iKoppenh oefer et al.' (2009) also used a Besangon model as pre- 
sented in Robin et al. (2003) to derive a brightness distribution of 
stars in the target field and performed Monte-Carlo simulations 
to simulate the occuiTence and detections of transits. These stud- 
ies include detailed observational constraints such as observing 
schedule, weather conditions, and exposure time and issues of 
data reduction, e.g. red noise and the impact of the instrument's 
point spread function. 

In our study, we rely on the extensive data reservoir of the 
Tycho catalog instead of assuming a stellar distribution or a 
Galactic model. We first estimate the number of expected ex- 
oplanet transit events as a projection on the complete celestial 
plane. We refer to recent results of transit surveys such as sta- 
tistical, empirical relationships between stellar properties and 
planetary formation rates. We then use basic characteristics of 
current low-budget but high-efficiency transit programs (BEST, 
XO, SuperWASP, and HATNet), regardless of observational con- 
straints mentioned above, and a simple model to test putative 
transits with the given instruments. With this procedure, we yield 
sky maps, which display the number of expected exoplanet tran- 
sit detections for the given surveys, i.e. the transit sky as it is 
seen through the eyeglasses of the surveys. 

The Tycho catalog comprises observations of roughly 1 mil- 
lion stars taken with the Hipparcos satellite between 1989 and 



^Extrasolar Planets Encyclopedia (EPE): www.exoplanet.eu. Four 
of these 62 announced transiting planets have no published position. 



1993 (lESAlll997trHoeglll997h . During the survey, roughly 100 
observations were taken per object. From the derived astromet- 
ric and photometric parameters, we use the right ascension (a), 
declination (6), the color index (B-V), the apparent visible mag- 
nitude my, and the stellar distance d that have been calculated 
from the measured parallax. The catalog is almost complete for 
the magnitude limit my < 11. 5", but we also find some fainter 
stars in the list. 



2. Data analysis 

The basis of our analysis is a segmentation of the celestial plane 
into a mosaic made up of multiple virtual fields of view (FsOV). 
In a first approach, we subdivide the celestial plane into a set of 
181 X 361 = 65 341 fields. Most of the current surveys do 
not use telescopes, which typically have small FsOV, but lenses 
with FsOV of typically 8° x 8°. Thus, we apply this extension of 
8° X 8° and a stepsize of A(5 = 1° = Aff, with an overlap of 7° 
between adjacent fields, for our automatic scanning in order to 
cover the complete sky. We chose the smallest possible step size 
in order to yield the highest possible resolution and the finest 
screening, despite the high redundancy due to the large overlap. 
A smaller step size than 1 ° was not convenient due to limitations 
of computational time. An Aitoff projection is used to fold the 
celestial sphere onto a 2D sheet. 

2. 1 . Derivation of tiie stellar parameters 

One key parameter for all of the further steps is the effective 
temperature Teif of the stars in our sample. This parameter is not 
given in the Tycho catalog but we may use the stellar color index 
(B - y) to deduce Tefj by 



r -f = 2Q[14.551-(B-y)]/3.684 j^ 



(1) 



which is valid for m ain-sequenc e (MS) stars with T^s < 9 100 K 
as late as type M8 (lReedlll998l) . Although we apply this equa- 
tion to each object in the catalog, of which a significant frac- 
tion might exceed Teff = 9 lOOK, this will not yield a serious 
challenge since we will dismiss these spurious candidates below. 
From the object's distance to Earth d and the visible magnitude 
my, we derive the absolute visible magnitude My via 



My 



OTv -5'"log 



lOpc 



(2) 



where we neglected effects of stellar extinction. In the next step, 
we compute the stellar radius in solar units via 



Re 



5 770K 



I eff 



2q(4.83-Mv)/2.5 



1/2 



and the stellar mass by 



(3) 



(4) 



where ctsb is the Stefan-Boltzmann constant. The coefficient /? 
in the relation L cc M^ depends on the stellar mass. We use 
the values and ma ss regimes that were empirically fo und by 
I Cester et all (Il983h . which are listed in Table [U (see also lSmithI 
;1983|k 

We deduce the stellar metallicity [Fe/H],^ from the star's ef- 
fective temperature Tetj and its color index {B - V) by 
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Table 1. Empirica l values for P in the mass-luminosity relation 
Eq. (IHi as given in lCester et"al1(ll983h . 





Stellar Mass Regime 


3.05 ±0.14 


< 0.5 Mo 


4.76 ±0.01 


O.6M0 <M* < 1.5 Mo 


3.68 ± 0.05 


1.5Mo <M* 



[Fe/H], = 



411 I K 



8423 + 4736(B- y)- 1 106(B- 



Vf) 
(5) 

as given in ISantos et al.l (l2004l) . This relation, however, is 
only vahd for stars with 0.51 < {B - V) < 1.33, 
4495K < Teff < 6339K, -0.7 < [Fe/H]* < 0.43, and 
\og{g) > 4. We reject those stars from the sample that do not 
comply with all these boundary conditions. On the one hand we 
cleanse our sample of non-MS stars, on the other hand the sam- 
ple is reduced seriously. While our original reservoir, our 'master 
sample', consists of 1 031 992 stars from the Tycho catalog, all 
the restrictions mentioned above diminish our sample to 392 000 
objects, corresponding to roughly 38 %. 

3. Transit occurrence and transit detection 

3.1. Transit occurrence 

Now that we derived the fundamental stellar parameters, we may 
turn towards the statistical aspects of planetary occurrence, ge- 
ometric transit probability and transit detection. We start with 
the probability for a certain star of the Tycho catalog, say the 
;* star, to host an exopla net. For F, G, and K d warfs with 
-0.5 < [Fe/H], < 0.5. iFischer & Valentil (l2005l) found the 
empirical relationship 



PBplanet,; = 0.03 ■ IQ^^^^/^^^ 



(6) 



for a set of 850 stars with an analysis of Doppler measurements 
sufficient to detect exoplanets with radial velocity (RV) semi- 
amp litudes K > 30ms ' and orbital periods shorter than 4yr 
(see lMarcv et al.l (l2005h and lWvatt eTal] (|20Q7l) for a discussion 
of the origin of this formula and its implications for planet for- 
mation). These periods are no boundary conditions for our sim- 
ulations since we are only interested in surveys with observing 
periods of < 50 d. The additional constraint on the metallicity 
does not reduce our diminished sample of 392 000 stars since 
there is no star with -0.7 < [Fe/H], < - 0.5 in the Tycho cat- 
alog. Similar to the correlation we use, lUdrv & SantosI (l2007h 
found a metalhcity distribution of exoplanet host stars equiv- 
alent to p3planet,i = 0.044 ■ 102"4[Fe/H], HoWCVCr, this fit 

was restricted to stars with [Fe/ HU > since the y suspect 
two regimes of planet f ormat ion. ISozzetti et al. j2009l) extended 

' (120051) 



the uniform sample of iFisch er & Va lentil (120051) and found the 
power-law jJapianet./ = 1-3 ■ 10- 1^"'")' +C,C e (0,0.5), to yield 
the best data fit. These recent studies also suggest that there ex- 
ists a previously unrecognized tail in the planet-metallicity dis- 
tribution for [Fe/H] ,^ < 0. Taking Eq. (|6]l we thus rather under- 
estimate the true occurrence of exoplanets around the stars from 
the Tycho catalog. The metallicity bias of surveys using the RV 
method for the detection of exoplanets is supposed to cancel out 
the b ias of transit surveys ( iGould et al.l 120061: iBeattv & Gaudil 
I2008h . 



In the next step, we analyze the probability of the putative 
exoplanet to actually show a transit. Considering arbitrary incli- 
nations of the orbital plane with respect to the observe r 's line 
of sight and including Kepler's third law. iGilUland et all (l2000h 
found the geometric transit probability to be 



^geo./ 



23.8 



Mo 



-1/3 



Rjid) 



-2/3 



(7) 



where P is the orbital period. A more elaborate expression - 
including eccentricity, planetary radius, the argument of perias- 
tron and the semi-major axis in stead of the orbital period - is 
given bv lSeagroves et all (l2003h . Note that pgeo./ in Eq. (|7]i does 
not explicitly but implicitly depend on the semi-major axis a via 
P - P(fl)! The probability for an exoplanetary transit to occur 
around the star is then given by 



— §^3 planet,/ ' Pgeo,/ 



(8) 



where P is the remaining free parameter, all the other parameters 
are inferred from the Tycho data. Since we are heading for the 
expectation value, i.e. the number of expected transits in a cer- 
tain field of view (FOV), we need a probability density for the 
distribution of the orbital periods of extrasolar planets. On the 
basis of the 233 e xoplanets listed in the EPF on July 6* 2007, 
iJiang et al.l (l 2007h used a power-law fit 6(P) = C{k) ■ (P/d)"*, 
with C{k) as the normalization function, and the boundary con- 
dition for the probability density dP 6(P) = 1 to get 



6{P)^ 



l-k 



B 



l-k 



\\-k 



0"' 



(9) 



with A = 1.21 1909d and B - 45 1 7.4 d as the lower and upper 
limits for the period distribution and k - 0.9277. This function 
is subject to severe selection effects and bases on data obtained 
from a variety of sur veys and instrumen ts. It overestimates short- 
period planets since I Jiang et al.l (|2007) included transiting plan- 
ets and the associated selection effects. While the function pre- 
sumably does not mirror the true distribution of orbital periods 
of exoplanets, it is correlated to the period distribution to which 
current instruments are sensitive, in addition to geometric selec- 
tion effects as given by Eq. (Q. 

We now segment the celestial plane into a mosaic made up 
of multiple virtual FsOV, as described at the beginning of Sect. 
|2] to calculate the number of expected transits in that field. In 
Sect. l3.2l we will attribute the FOV of the respective instrument 
to that mosaic and we will also consider the CCD resolution. The 
number of stars comprised by a certain FOV is n. The number of 
expected transits around the star in that field, A^,, with periods 
between Pi and P2 is then given by 



A^, = dP d(P) I 



(10) 



M, \ ^'^ I Ri\ 1 -k 



y 1 (pUi-k pl/3-nH2/3 

1/3 -A: '2 1^ \A<p,<P2<B , 

and the number of expected transits in the whole FOV is 

n 

1=1 



(11) 
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1 2 3 4 5 7 8 91011 
Expected No. of Transit Events 




-200 -100 100 

an 



Fig. 1. Sky map of the expected number of exoplanet tran- 
sit events, A^, with orbital periods between Pi = 1.5 d and 
P2 = 50 d on the basis of 392000 objects from the Tycho cat- 
alog. The published positions of 58 transiting planets from the 
EPE as of September 1'*' 2009 are indicated with symbols: 6 de- 
tections from the space-based CoRoT mission are labeled with 
triangles, 52 ground-based detections marked with squares. The 
axes only refer to the celestial equator and meridian. 



We emphasize that this is not yet the number of expected transit 
detections within a certain FOV (see Sect. [Hi but the number of 
expected transits to occur within it. 

A graphical interpretation of this analysis is presented in Fig. 
[1] where we show a sky map of the expected number of exo- 
planet transits around MS stars with my < 11 .5" for orbital 
periods between Pi = 1.5dandP2 = 50 d. This map bases on 
several empirical relationships and on substantial observational 
bias towards close-in Jupiter-like planets, but nevertheless it rep- 
resents the transit distribution to which current instrumentation 
has access to. The pronounced bright regions at the upper left 
and the lower right are the anti-center and the center of the Milky 
Way, respectively. The absolute values of 0.5 < N < 5 for the 
most of the sky are very well in line w ith the experiences from 
wide-field surveys using a 6° x 6° field. iMandushev et al.l ( l2005h 
stated 5 to 20 or more exoplanet transit candidates, depending on 
Galactic latitude, and a ratio of =s 25 : 1 between candidates 
and confirmed planets, which is equivalent to 0.2 < < 1. 
Our values are a little higher, probably due to the slightly larger 
FOV of 8° X 8° used in Fig. [T]and due to the eff'ect of blends 
and unresolved binaries (see discussion in Sect. [Sj. 

In the left panel of Fig. |2] we show the distribution of ex- 
pected transits from our simulation as a function of the host 
stars' magnitudes compared to the distribution of the observed 
transiting exoplanets. The scales for both distributions differ 
about an order of magnitude, which is reasonable since only a 
fraction of actual transits is observed as yet. For my < 8™, only 
HD209458b and HD 189733b are currently know to show tran- 
sits whereas we predict 20 of such transits with periods between 
1.5 d and 50 d to occur in total. We also find that the number of 
detected transiting planets does not follow the shape of the sim- 
ulated distribution for my > 9™. This is certainly induced by 
a lack of instruments with sufficient sensitivity towards higher 
apparent magnitudes, the much larger reservoir of fainter stars 
that has not yet been subject to continuous monitoring, and the 
higher demands on transit detection pipelines. 



Our transit map allows us to constrain convenient locations 
for future ground-based surveys. A criterion for such a loca- 
tion is the number of transit events that can be observed from 
a given spot at latitude / on Earth. To yield an estimate, we in- 
tegrate over that part of the celestial plane that is accessable 
from a telescope situated at I. We restrict this observable fan to 
/ - 60° < 6 < I + 60°, implying that stars with eleva- 
tions > 30° above the horizon are observable. The number of 
the transit events with my < 1 1.5™ that is observable at a cer- 
tain latitude on Earth is shown in the right panel of Fig.|2l This 
distribution resembles a triangle with its maximum almost ex- 
actly at the equator. Its smoothness is caused by the wide angle 
of 120° that flattens all the fine structures that can be seen in Fig. 

m 

3.2. Transit detection 

So far, we have computed the sky and magnitude distributions 
of expected exoplanet transits with orbital periods between 1.5 d 
and 50 d, based on the stellar parameters from the Tycho data 
and empirical relations. In order to estimate if a possible tran- 
sit can actually be observed, one also has to consider technical 
issues of a certain telescope as well as the efficiency and the se- 
lection effects of the data reduction pipelines. The treatment of 
the pipeline will not be subject of our further analysis. The rel- 
evant aspects for our concern are the pixel size of the CCD, its 
FOV, the my range of the CCD-telescope combination, and the 
declination fan that is covered by the telescope. 

To detect a transit, one must be able to distinguish the pe- 
riodic transit pattern within a light curve from the noise in the 
data. Since the depth of the transit curve is proportional to the 
ratio Ap/Ai,, where Ap and A* are the sky-projected areas of the 
planet and the star, respectively, and Rp/Ri, = yfApjA^, with Rp 
as the planetary radius, the detection probability for a certain in- 
strument is also restricted to a certain regime of planetary radii. 
Assuming that the transit depth is about 1 %, the planetary ra- 
dius would have to be larger than ^ /J*/10. We do not include 
an elaborate treatment of signal-to-noise in our considerations 
(see Aigrain & Pont 2007':). Since our focus is on MS stars and 
our assumptions for planetary occurrence are based on those of 
Hot Jupiters, our argumentation automatically leads to planetary 
transits of exoplanets close to x P^./10. 

We also do not consider observational aspects, such as in- 
tegration time and an observer on a rotating Earth with ob- 
servation windows a nd a finite amount of observing time (see 
[Fleming et al.l l20()8l for a review of these and other observa- 
tional aspects). Instead, we focus on the technical characteristics 
of four well-established transit surveys and calculate the celestial 
distribution of expected exoplanet transit detections in principle 
by using one of these instruments. The impact of limited observ- 
ing time is degraded to insignificance because the span of orbital 
periods we consider in Eq. (fTOl i reaches only up to P2 = 50 d. 
After repeated observations of the same field, such a transiting 
companion would be detected after < 3 yr, which is the typical 
duty cycle of current surveys. 

Our computations are compared for four surveys: BEST, 
XO, SuperWASP, and HATNet. This sample comprises the three 
most fruitful surveys in terms of first planet detections and BEST 
- a search program that used a telescope instead of lenses. 
While observations with BEST have been ceased without any 
confirmed transit detection, XO has announced detections and 
SuperWASP and HATNet belong the most fruitful surveys to 
date. An overview of the relevant observational and technical 
properties of these surveys is given in Table |2] For each survey. 
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Fig. 2. Sky-integrated number of transits per magnitude (left panel) and as a function of latitude (right panel). Left: While the green 
line represents our simulations, the rosy bars show the number of transiting planets per magnitude bin discovered so far. Note the 
different scales at the left and right ordinates! Right: The triangle represents the expected number of transits that can be seen at 
elevations higher than 30° over the horizon at a given latitude on Earth. 



we first restrict the Tycho master sample to the respective mag- 
nitude range, yielding an ;7Zv -restricted sample. In the next step, 
we virtually observe the subsample with the fixed FOV of the 
survey telescope, successively grazing the whole sky with steps 
of 1° between adjacent fields. The FOV is composed of a num- 
ber of CCD pixels and each of these pixels contains a certain 
number of stars, whose combined photon fluxes merge into a 
count rate. Efficient transit finding has been proven to be pos- 
sible from the ground in crowded fields, where target objects 
are not resolved from neighbor stars. To decide whether a hy- 
pothetical transit around the /* star in the pixel would be de- 
tected, we simulate the effect of a transiting object that reduces 
the light flux contribution /, of the star on the combined flux 
Y!l h of the stars within a pixel. If the /* magnitude variation on 
the pixel-combined light is Amy,, > 0.01", which is a typical 
accuracy limit of current ground-based surveys, then we keep 
this star for further analysis of the transit detection as described 
in Sects. I2land l3. II otherwise it is rejected. The fluxes, however, 
are not listed in the Tycho catalog; instead, we can use the vis- 
ible magnitude wzy,/ of a star and calculate its relative flux ///o 
with respect to a reference object with flux /o at magnitude otv,o: 



fi 



(12) 



The magnitude variation can then be computed via 



0.01 < Amy,/ = - 2.5 ■ log 



Q.99- fi + Yj h 

k 



(13) 



2.5 ■ log 



0.99. f.^f^ 



Zfk_ 
fo 



Without loss of generality we chose my.o = 30™ as reference 
magnitude. 



4. Results 

We develop a procedure for the calculation of the number of 
expected transit events to occur around MS stars based on em- 
pirical relations between the stars and planets. This procedure is 
then applied to more than 1 million stars from the Tycho catalog 
to visualize the transit probability for all stars with my < 11 .5™ 
as a sky map. We also compute the celestial distribution of the 
number of expected transit detections for four different, well- 
established wide-field surveys. 

In Fig. [3] we present the number of expected transit detec- 
tions for the technical properties of BEST, XO, SuperWASP, 
and HATNet. As a general result from these maps, we find that 
the size of the FOV governs the detection efficiency of a cam- 
era. For the method applied here, the CCD resolution, i.e. the 
pixel size, has almost no impact since we neglect effects of noise, 
whereas in general the detection limits for transiti ng planets de- 
pend on the CCD resolution in terms of noise dKovacs et al.l 
120051: iTamuz et alj|2005l: iPont et al.ll2006l) . In Table |2] you see 
that the restriction of Amy > O.OI™ almost doesn't reduce 
the sample. A large FOV, collecting the light of relatively many 
stars, outweighs a lower CCD resolution - at least for the range 
of pixel sizes considered here. Even for the zones around and 
in the galactic center and anti-center where the stellar density 
increases drastically, the number of detectable transit events 
reaches its maximum. This was not foreseeable since blending, 
simulated by Eq. ( fT3b , could have reduced the efficiency of tran- 
sit detection within the crowded zones. 

The four survey sky maps portray the very distinct efficien- 
cies of the telescopes. The map of BEST reflects the stellar dis- 
tribution of the Tycho data best due to the relatively high reso- 
lution of the CCD. However, the small FOV leads to very few 
expected transit detections. BEST'S visible magnitude cut at the 
upper end is 14'" while the Tycho catalog is complete only up to 
11.5™. Thus, a significant contribution of stars inside this range 
is excluded in Fig.|3^. BEST also covers the smallest portion of 
the sky, compared to the other surveys. The XO project yields 
a much more promising sky map, owed to the larger FOV of 
the lenses. But due to the relatively large pixel size and an ad- 
verse magnitude cut of 9 < my, XO achieves lower densi- 
ties of expected detections than SuperWASP and HATNet. As 
for SuperWASP and HATNet, the difference in the magnitude 
cuts with respect to the Tycho catalog is negligible for XO but 
tends to result in an underestimation of the expected detections. 
That part of the SuperWASP map that is also masked by HATNet 
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looks very similar to the map of the latter one. While HATNet 
reaches slightly higher values for the expected number of detec- 
tions at most locations, the covered area of SuperWASP is sig- 
nificantly larger, which enhances its efficiency on the southern 
hemisphere. 

The total of expected transiting planets in the whole sky is 
3412 (see Fig. [TJ. By summing up all these candidates within 
an observational fan of 30° elevation above the horizon, we lo- 
calize the most convenient site on Earth to mount a telescope 
for transit observations (see right panel in Fig. |2]i: it is situated 
at geographical latitude / - -1°. Given that the rotation of the 
Earth allows a ground-based observer at the equator, where both 
hemispheres can be seen, to cover a larger celestial area than 
at the poles, where only one hemisphere is visible, this result 
could have been anticipated. Due to the non-symmetric distribu- 
tion of stars, however, the shape of the sky-integrated number of 
expected transits as a function of latitude is not obvious. Figure 
|2] shows that the function is almost symmetric with respect to 
the equator, with slightly more expected transits at the northern 
hemisphere. Furthermore, the number of expected transits to be 
observable at the equator is not twice its value at the poles, which 
is due to the inhomogeneous stellar distribution. In fact, an ob- 
server at the equator triples its number of expected transits with 
respect to a spot at the poles and can survey almost all of the 
3412 transiting objects. 

Based on the analysis of the magnitude distribution (left 
panel in Fig. [2]), we predict 20 planets with my < 8 to 
show transits with orbital periods between 1.5 d and 50 d, 
while two are currently known (HD209458b and HD 189733b). 
These objects have proven to be very fru itful for follow-up 
studies such as transmission spectro scopv ([Charbonneau et al 



2002; 



2008 



Wlal-Madiar et al. 2004; Knu tson et al.ll2007l: ISing et al 



Grillmair et al.n2008c iPont et alj|2008h easurements 
of the Rossiter-McLaughUn effect (lHoltlll893HRossiteij 1 19241: 
McLaughlin' 192?; Winn et al. '2006', '2007; Wolfet al.' 7007, 
Narita et al.. .2007; Cochran et al. 2008 ; Winn et aL 200 9). Our 
analysis suggests that a significant number of bright transiting 
planets is waiting to be discovered. We localize the most promis- 
ing spots for such detections. 



5. Discussion 

Our valu es for the XO project are much higher than those pro- 
vided by iBeattv & Gaudi, (2008), who also simulated the ex- 
pected exoplanet transit detections of XO. This is due to their 
much more elaborate inclusion of observational constraints such 
as observational cadence, i.e. hours of observing per night, me- 
teorologic conditions, exposure time, and their approach of mak- 
ing assumptions about stellar densities and the Galactic structure 
instead of using catalog-based data as we did. Given these differ- 
ences between their approach and ours, the results are not one- 
to-one comparable. While the study of iBeattv & Gaudil (l2008h 
definitely yields more realistic values for the expected number 
of transit detections considering all possible given conditions, 
we provide estimates for the celestial distribution of these detec- 
tions, neglecting observational aspects. 

In addition to the crucial respects that make up the effi- 
ciency of the projects, as presented in Table |2] SuperWASP and 
HATNet benefit from the combination of two observation sites 
and several cameras, while XO also takes advantage of twin 
lenses but a single location. Each survey uses a single camera 
type and both types have similar properties, as far as our study is 
concerned. The transit detection maps in Fig. [3] refer to a sin- 
gle camera of the respective survey. The alliance of multiple 



cameras and ffie diverse observing strategies arn ong the surveys 
dMcCuUough et al.ll2005t ICameron et al.ll2009l) bias the speed 
and efficiency of the mapping procedure. This contributes to the 
dominance of SuperWASP (18 detections, 14 of which have pub- 
lished positionsJJ over HATNet (13 detections, all of which have 
published positions)^, XO (5 detections, all of which have pub- 
lished positions)'', and BEST (no detection)^. 

It is inevitable that a significant fraction of unresolved bi- 
nary stars within the Tycho data blurs our results. The impact of 
unresolved binaries without physical interaction, which merely 
happen to be aligned along the line of sight, is significant only 
in the case of extreme crowding. As shown bv lGillon & MagainI 
(l2007h . the fraction of planets not detected because of blends is 
typically lower than 10%. The influence of unresolved physi- 
cal binaries will be higher. Based on the empirical peri od dis - 
tribution for binary s tars from Duquennov & Mavor ( 199lh . 
IBeattv & Gaudil (l2008l) estimate the fraction of transiting plan- 
ets that would be detected despite the presence of binary systems 
to be » 70%. Both the contribution of binary stars aligned by 
chance and physically interacting binaries result in an overesti- 
mation of our computations of » 40 %, which is of the same 
order as uncertaintie s arising from the e mpirical relationships 
we use. Moreover, as IWillemsetal](l2006l) have shown, the den- 
sity of eclipsing stellar binary systems increases dramatically to- 
wards the Galactic center. To control the fraction of false alarms, 
efficient data reduction pipe lines, and in particular data analysis 
algor ithms, are necessary dSchwarzenberg-Czernv & Beaulieul 
I2006h . 

Recent evidence for the existence of ultra -short period plan- 
ets around low-mass stars dSahu et al.ll2009l) . with orbital peri- 
ods < Id, suggests ffiat we underestimated the number of 
expected transits to occur, as presented in Sect. 13.11 The possi- 
ble underestimation of exoplanets occurring at [Fe/H]^, < 
also contributes to a higher number of transits and detections 
than we computed here. Together with the fact that the Tycho 
catalog is only complete to my < 11.5™, whereas the surveys 
considered here are sensitive to slightly fainter stars (see Table 
|2]i, these trends towards higher numbers of expected transit de- 
tections might outweigh the opposite effect of unresolved binary 
stars. 

A radical refinement of both our maps for transits occurrence 
and detections will be available within the next few years, once 
the 'Panoramic Survey Telescope a nd Rapid Response System' 
(Pan-STARRS) (iKaiser et al.ll2002l) will run to its full extent. 
Imaging roughly 6000 square degrees every night with a sensi- 
tivity down to my x 24, this survey will not only drastically in- 
crease the number of cataloged stars - thus enhance our knowl- 
edge of the localization of putative e xoplanetary transits - but 
could potentially detect transits itself dPupuv & Liull2009h . The 
Pan-STARRS catalog will provide the ideal sky map, on top of 
which an analysis presented in this paper can be repeated for 
any ground-based survey with the aim of localizing the most ap- 
propriate transit spots on the celestial plane. The bottleneck for 
the verification of transiting planets, however, is not the local- 
ization of the most promising spots but the selection of follow- 
up targets accessible with spectroscopic instruments. The ad- 
vance to fainter and fainter objects thus won't necessarily lead 
to more transit confirmations. Upcoming spectro graphs, such as 
the E SPRESSO@VLT and the CODEX@E-ELT dPepe & LovisI 
l2008l) . can be used to confirm transits around fainter objects. 
These next-generation spectrographs that will reveal Doppler 
fluctuations on the order of cm s ' will also enhance our knowl- 
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Fig. 3. Sky maps with expected number of transit detections for BEST (a), XO (b), SuperWASP (c), and HATNet (d). Note the 
different scales of the color code! The published detections of the surveys, taken from the EPE as of September P' 2009, are 
indicated with symbols. 



edge about Hot Neptunes and S uper-Earths, whic h the recently 
discovered transits of GJ 436 b ( iButler et al.ll2004h. HA T-P- 1 1 b 
dBakos etal.ll2009h. and CoR oT-7b dLeger et al.ll2009l) and re- 



Lovis et al] (l2009t) predict to be numerous 



suits from 

Further improvement of our strategy will emerge from the 
findings of more exoplanets around MS stars and from the us- 
age of public data reservoirs like the NASA Star and Exoplanet 
DatabaseQ, making assumptions about the metallicity distribu- 
tion of planet host stars and the orbital period distribution of 
exoplanets more robust. 
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Table 2. Instrumental properties of the treated surveys. 



Survey 


S RANGE 


FOV 


CCD Pixel Size 


m\i Range 


6 & mv 


S, my & Amy 


S, my, Amy & MS 










[Mag.] 


Limited Sample 


Limited Sample 


Limited Sample 


BEST 


-16° < 5 < 90° 


3.1° 


5.5 "/mm 


8 < OTv < 14 


546 382 (52.94%) 


516 524 (50.05%) 


222 854 (21.59%) 


XO 


-39° < a < 90° 


7.2° 


25.4 "/mm 


9 < mv < 12 


620 477 (60.12%) 


597 842 (57.93%) 


263 213 (25.51%) 


SuperWASP 


-44° < 5 < 90° 


7.8° 


13.8"/mm 


1 <m-^ < 12 


745 227 (72.21 %) 


703 707 (68.19%) 


311404 (30.18%) 


HATNet 


-28° < 5 < 90° 


8.3° 


14.0 "/mm 


7 < OTv < 12 


721 473 (69.91 %) 


686927 (66.56%) 


283 350 (27.46%) 



In the last three columns we list the reduced Tycho master sample of 1 031 992 stars after we applied the subsequent boundary conditions: the 
survey's sky-coverage (d range), its my limitation, magnitude variation Amy > 0.01" for the transit of a Jupiter-sized object around the star in 
a pixel, and the boundary conditions for MS stars, for which the empirical relationships hold (see Sect. 12. U . In braces we indicate the portion of 
the Tycho master sample. 



